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ABSTRACT 


An experimental investigation dealing with the fracture 
behavior of unidirectional boron/aluminum composite lami- 
nates is conducted in order to verify the results of mathe- 
matical models developed by the author and associates. 

These models predict the stresses and displacements of 
fibers and the amount of damage growth in a center-notched 
lamina as a function of the applied remote stress and the 
matrix and fiber material properties. The damage may take 
the form of longitudinal yielding and splitting in the 
matrix as well as stable transverse damage consisting of 
broken fibers and matrix yielding ahead of the notch. A 
brittle lacquer coating is used to detect the yielding in 
the matrix while X-ray techniques are used to determine the 
number of broken fibers in the laminate. 

The notched strengths and the amounts of damage found 
in the experimental specimens agree well with those pre- 
dicted by the mathematical model. Results of the tests are 
presented in the form cf graphs and photographs. It is 
shown that the amount of damage and the fiber displacements 
do not depend strongly on the number of plies in the lami- 
nate for a given notch width. By heat-treating certain 
laminates to increase the yield stress of the aluminum 
matrix, the effect of different matrix properties on the 
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fracture behavior is investigated. It is demonstrated that 
the stronger matrix will actually weaken the notched lami- 
nate by decreasing the amount of matrix damage and thereby 
making the laminate more notch sensitive. 
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CHAPTER I 


INTRODUCTION 

The desire for lighter and stronger structural materi- 
als has lead to a rapidly growing use of fibrous composite 
materials in many modern applications. By combining the 
properties of more than one material having different 
strengths and weights; a structural member that is lighter 
and stronger than the isotropic material it replaces can be 
produced. An added benefit often gained is that due to this 
multiphase construction composites are frequently more 
insensitive to damage and damage growth than isotropic mate- 
rials. 

Many different types of composite materials can be 
found in use today such as short fiber composites, woven 
cloth composites, and composites with randomly placed 
fibers. The list of matrix materials used to support these 
fibers is at least as long as the list of different fibers. 
In advanced structures such as airframes, however, most 
structural members are made of multiple layers, each con- 
sisting of continuous, parallel fibers embedded in a matrix. 
These layers, each of which is referred to as a lamina, are 
bonded together to form the structural member, typically 
called a laminate. 
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By orienting these continuous fibers in specific direc- 
tions the designer can adjust the material properties of the 
laminate to suit the loading to be encountered in use. Some 
composites having the predominate loading in a given direc- 
tion, such as landing struts and wing panels, are often con- 
structed of all unidirectional composite laminates. Others, 
which encounter a more randomly oriented load spectrum, are 
made by constructing a laminate from a series of properly 
oriented unidirectional lamina. 

In addition to the fiber orientation, advanced compos- 
ite materials are further grouped into two broad categories 
based on the type of matrix material used, these being plas- 
tic matrix composites and metal matrix composites. Unidi- 
rectional composites which have epoxy (plastic) matrices 
characteristically form largo splits in the matrix between 
fibers near a defect such as a notch. Ductile metal matrix 
composites, on the other hand, exhibit matrix yielding 
between fibers near a notch instead of this splitting. It 
is found to be possible to have stable notch extension and 
broken fibers in composites which use ductile Ktetals such as 
aluminum for the matrix material, presumably because of this 
large scale yielding. This is one question to be considered 
in the present study. 

For both unidirectional and multidirectional composites 
of the type discussed above, the fundamental component is 
the unidirectional lamina. It is thus essential that the 
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behavior of the lamina be clearly understood before the 
fracture behavior of the multiply laminate can be described. 
Many of the mathematical models which have been developed 
consider this fundamental lamina and the present experimen- 
tal study will attempt to relate the predictions of some of 
these models to actual laboratory measurements. 

Considerable analytical work has been done in modeling 
the unidirectional lamina with many of the recent mathemati- 
cal models developed at Clemson University by the authors 
[1,2). These models, based on earlier work by Hedgepeth [3] 
in 1961, assume that the parallel fibers in the lamina carry 
all of the axial load and the matrix material transmits load 
to adjacent fibers by shear. An influence function techni- 
que was used by Hedgepeth to solve the lamina equilibrium 
equations. In later work, Hedgepeth and Van Dyke [4] were 
able to include matrix yielding and splitting in their 
model, although they were restricted by the solution techni- 
que to only one broken fiber in the lamina. Later models, 
such as the model of Eweben [5], extended Hedgepeth's model 
to include an arbitrary number of broken fibers while still 
allowing for matrix damage. Zweben, however, simplified the 
model considerably by grouping the fibers in three regions, 
one including all of the broken fibers, one having unbroken 
fibers which are influenced by the notch, and one group of 
unbroken fibers which are not influenced by the defect. 
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A limited amount of experimental results was available 
at the time of publication of [1], with perhaps the most 
complete being that of Awerbuch and Hahn [6] using narrow 
unidirectional boron/aluminum panels, These center-notched 
panels were loaded to failure so that their notched 
strengths could be determined. Extensive monitoring of the 
crack opening displacement (COD) was performed during the 
tests, along with some observations concerning the stable 
notch extension and matrix damage. This laminate damage was 
monitored through the use of radiographs and photomicro- 
graphs. The COD, defined as twice the vertical displacement 
of the center fiber, was monitored through the use of an 
interferometric displacement gage. 

Awerbuch and Hahn were able to predict the load-COD 
curves using a simplified analytical model which could 
account for some crack-tip yielding in the matrix. Results 
for a narrow range of initial notch widths were presented 
along with some interpretation of the relation between the 
observed failure modes and deformation characteristics. All 
of the tests performed by Awerbuch and Hahn, however, were 
carried out using narrow, 25.4 mm wide specimens. A finite 
width correction factor was used in the study to allow the 
predicted stresses for the infinite panel to be compared 
with the narrow specimens. These investigators assumed that 
the finite width correction factors available for isotropic 
plates would be applicable to their fibrous composite 
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specimens. But even with the finite width correction fac- 
tors, these investigators were limited to notch widths of 
12.7 mm (roughly 71 broken fibers) or less because of the 
narrow specimen width. All of the, specimens used in the 
experiments were eight-ply; thus nothing was said concerning 
the effect of laminate thickness on the damage growth and 
failure mechanisms. 

Grandemange and Street [7] performed tests on double 
edge-notched boron/aluminum panels in an early attempt to 
investigate fracture mechanisms in fibrous composite materi- 
als. Their investigation included unidirectional laminates 
having a reasonably wide range of thicknesses, but a very 
narrow range on notch widths. The suggestion was made by 
these authors that the thickness of a specimen did affect 
the notched strength of an edge-notched panel. However, the 
strength of a panel having a given notch width varied only 
10 percent over the entire range of thicknesses tested. It 
is clear that a much wider range of notch widths must be 
tested before the thickness effect can be thoroughly 
explained. 

Additional experimental work has been performed more 
recently by Reedy [8] also using boron/aluminum laminates, 
where very thick (28 ply) specimens were center-notched and 
loaded to failure. Reedy suggested in his paper that the 
thick specimens he used for his experiments should have been 
more notch sensitive than monolayer panels which he modeled 
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analytically. Yet in the results published by Reedy, the 
monolayer model predictions appeared to agree with the expe- 
rimental results for the thick specimens reasonably well. 
However, as in the work of Awerbuch and Hahn [6], only one 
thickness was used in the experiments, and the question of 
thickness effect on fracture behavior was not addressed. 

Reedy also noted the large amounts of longitudinal 
matrix yielding predicted by the analytical model, and 
attempted to measure this yielding in the test laminate at 
various load levels using strain gage rosettes mounted at 
various points on the panel. The yield region predicted 
analytically appeared to agree fairly well with the strains 
measured in the laboratory, but since the strain gages used 
were 3 mm in width and thus covered nearly 20 fibers, there 
is some question as to the accuracy of these readings. The 
gages were clearly connected to some fibers that were broken 
and some that were not, thus the output of the gages needs a 
more careful analysis than is available from standard 
stress-strain relations, 

Other attempts have been made to describe the damage 
growth in fibrous composite panels during fracture. Rollins 
[9] and Grenis and Levitt [10] attempted to relate the fiber 
breaks during a loading cycle to the acoustic emissions 
detected by transducers mounted on the specimens. It was 
found that the acoustic emission associated with a fiber 
break could be detected since these emissions were quite 
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energetic compared to those caused by matrix yielding and 
splitting. The points at which the various damage modes 
would occur in the loading sequence and the damage necessary 
to break a fiber were, however, never adequately explained. 

The need for additional experimental results became 
apparent during the final stages of the analytical modeling 
by Dharani, Jones, and Goree [11]. It was felt that a clear 
and logical description of the damage growth and fracture 
behavior in unidirectional laminates was provided by the 
micromechanical model. No detailed experimental studies 
describing this damage growth and fracture behavior were 
available, however, and much of the motivation for the work 
to be presented in this study was the need to observe very 
carefully the damage growth in a unidirectional composite 
and to compare the results with those predicted by the anal- 
ysis from [1] and [2], 

The analysis given in [1] and [2] is valid for both 
brittle (plastic) matrix and ductile (metal) matrix compos- 
ites with the fracture behavior of a brittle matrix laminate 
such as graphite/epoxy being relatively simple to describe. 
At very low loads, the brittle matrix splits at the ends of 
the notch with the splits soon reaching the ends of the 
specimen for little increase in load. The specimen then 
fails in net-section tension with essentially no stress con- 
centration effect other than the loss in area due to the 
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On the other hand the ductile matrix composite such as 
boron/aluminum exhibits large longitudinal yielding and 
transverse broken fibers and is thus much more complicated 
and difficult to model accurately. It was thus proposed 
that an experimental investigation be conducted using unidi- 
rectional boron/aluminum composite panels. It was known 
from past experience that finished boron/aluminum panels 
were available that possessed very uniform fiber size and 
spacing, and that the high quality of these finished speci- 
mens would make the tests much more valid. Very detailed 
records would be kept concerning the amount of matrix yield- 
ing or splitting as well as the number of broken fibers 
present in the laminate at various points in the loading. 

The results of these tests could then be compared to the 
damage zones predicted by the mathematical model. 

Unlike earlier tests, this proposed investigation would 
incorporate specimens of various widths and thicknesses 
(numbers of plies). Through the use of wider specimens, a 
much broader range of initial notch widths could be investi- 
gated than that of Awerbuch and Hahn [6]. These wider 
notchs could be cut in much wider specimens, reducing the 
need for the use of the finite width correction factor used 
in [6]. If the notchs in the proposed tests grew such that 
a finite width correction factor was needed, a much more 
accurate solution, recently been developed by Dharani [12] 
that is appropriate for the unidirectional composite 
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laminate, could be used. Also, the question of thickness 
effect on damage growth proposed in [7] and [8] could con- 
ceivably be answered through the use of specimens having 
numbers of plies ranging from one to eight. 

As in some earlier investigations, the number of broken 
fibers present in the laminate at various load levels could 
be determined using X-rays to penetrate the specimens. The 
X-ray records would be obtained using a portable X-ray 
source which would allow radiography while the specimen is 
in the testing machine under load, rather than unloading and 
placing the specimen in an X-ray cabinet as was done in [6]. 
It was hoped that photoelastic coatings could be used to 
monitor the amount of matrix yielding present in the lami- 
nates, since this method had been used by others such as 
Yeow, Morris, and Brinson [13]. But since the thinnest 
available coatings were thicker than many of the proposed 
specimens, it was clear after some investigation that the 
photoelastic coating would not be able to detect the yield- 
ing between fibers with the desired accuracy. 

It was thus decided that a brittle lacquer surface 
coating would be the best method by which to measure the 
yielding present in the matrix between the fibers. It had 
already been shown by Chaturvedi and Agarwal [14] that a 
brittle coating would crack at a given threshold strain in a 
direction perpendicular to the maximum principal strain even 
on a unidirectional composite material. Because the brittle 
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coating can be applied in a very thin layer it is sensitive 
to local effects which occur over dimensions on the order of 
the fiber spacing. Thus the location and orientation of the 
cracks on the specimens in the proposed tests should enable 
the investigators to determine the extent of yielding pres- 
ent at selected load levels. 

One aspect of damage growth in composites that has not 
been investigated thoroughly is the effect of different 
matrix properties on the damage growth and fracture behav- 
ior. Only some very recent work by Reedy [15,16] has 
attempted to determine the result of using a stronger or 
weaker aluminum matrix material to hold the fibers. This 
work was again limited in scope, using only one notch width 
and narrow, 25.4? mm wide specimens. It was thus suggested 
that some of the proposed specimens have an aluminum matrix 
that had been heat-treated to a T6 condition instead of the 
fully annealed 6061 aluminum matrix in the remainder of the 
panels. This stronger matrix should alter the fracture 
behavior, since it is through the matrix that the load is 
transferred from the broken fibers to the remaining unbroken 
ones. Furthermore, unlike the work in [15] and [16], a wide 
range of notch widths and panel thicknesses could be tested 
to provide much more information about the effect of matrix 
properties . 

A series of tests using a wide range of specimens hav- 
ing different notch widths, overall widths, numbers of 
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plies, and matrix properties was thus proposed. A complete 
description of the materials and test procedure used in the 
tests is described in the next section. 


CHAPTER II 


EXPERIMENTAL PROCEDURE 
Materials 

The specimens used in this investigation were all uni- 
directional boron/aluminum composite panels of various 
widths and thicknesses. The panels, manufactured by Amer- 
com, Inc. in Chatsworth, California, consisted of 0.142 mm 
(5.6 mil) diameter boron fibers in a 6061 aluminum matrix. 
Some of the panels were heat-treated after the initial manu- 
facture such that the aluminum matrix would possess the 
properties of standard 6061-T6 aluminum instead of an ann- 
ealed 6061-0 aluminum matrix. Equal numbers of one-, two-, 
four-, and eight-ply specimens were fabricated, with a 
thickness of approximately approximately 0.178 mm per ply. 

A complete listing of the test specimens used in this study 
can be found in Table (1). It should be noted that a set of 
specimens is defined as one single-ply specimen, one two-ply 
specimen, one four-ply specimen, and one eight-ply specimen, 
all having the same notch size and laminate width. Slightly 
thicker aluminum foil was used on the single-ply panels in 
order to increase out-of-plane stiffness. Thus the fiber 
volume fraction for the single-ply panels was around 35 per- 
cent, while the volume fraction for all other specimens was 
very close to 50 percent. 
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Table 1. List of specimens used in experimental 
investigation. 

Panel Width (W) Notch Width (2a) Matrix Temper* 


25. 4 

mm 

1.59 

mm 



1 set 6061-0 


25.4 

mm 

3.18 

mm 

i 

set 

6061-0, 1 set 

6061-T6 

25.4 

mm 

4.76 

mm 

i 

set 

6061-0, 1 set 

6061-T6 

25.4 

mm 

6.35 

mm 



1 set 6061-0 


50,8 

mm 

1.59 

mm 



1 set 6061-0 


50.8 

mm 

6.35 

mm 

i 

set 

6061-0, 1 set 

6061-T6 

50,8 

mm 

9.52 

mm 

i 

set 

6061-0, 1 set 

6061-T6 

50.8 

mm 

15.9 

nun 



1 set 6061-0 


73.1 

mm 

3.18 

mm 



1 set 6061-0 


73.1 

mm 

12.7 

mm 

i 

set 

6061-0, 1 set 

6061-T6 

73.1 

mm 

15.9 

mm 

i 

set 

6061-0, 1 set 

6061-T6 

73.1 

mm 

19.1 

mm 



1 set 6061-0 


101. 

6 mm 

15.9 

mm 



1 set 6061-0 


101. 

6 mm 

19.1 

mm 

i 

set 

6061-0, 1 set 

6061-T6 

101. 

6 mm 

22.2 

mm 



1 set 6061-0 


101. 

6 mm 

25.4 

mm 

i 

set 

6061-0, 1 set 

6061-T6 


* 

One set equals one single-ply panel, one two-ply 
panel, one four-ply panel, and one e:Lght-ply panel. 
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Experimental specimens having wiaths of 25.4, 

50.8,73.1, and 100.6 mm were cut from large panels obtained 
from the manufacturer. These cuts were parallel to the 
fibers and were made with a diamond saw. The edges were 
later ground to ensure minimum surface damage due to cut- 
ting. All of these experimental specimens were approxi- 
mately 30 cm in length. Center notches of various widths 
were then introduced by E.D.M. (Electrostatic Discharge 
Machining) . Normally the width of the notch did not exceed 
one-quarter of the overall width of the panel. 

Finally, the experimental specimens were carefully 
cleaned and polished to eliminate further the possibility of 
surface flaws that might affect the testing. Two strain 
gages were mounted on each specimen at points well below the 
notch so that the remote strain could be measured and 
recorded. The strain gages also enabled the experimenters 
to determine if the panel was indeed being subjected to a 
uniform remote stress. The last step in the preparation was 
to epoxy aluminum doublers to each end of each panel so that 
the testing machine grips would not damage the ends of the 
specimens. A representation of a typical test specimen is 
shown in Figure (1). 



Figure 1. Typical Test Specimen. 
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Radi.oqraphlc Technique 

The radiographs used in these experiments were taken by 
exposing the specimens to low level X-rays for a specified 
amount of time. The source of these X-rays was a Model MTK 
140Be X-ray machine manufactured by the Philips Company of 
West Germany. This is a portable unit and has the signifi- 
cant advantage that the radiographs can be taken while the 
specimens are in the testing machine and under load. 

Much tlmfr and experimentation was involved in determin- 
ing the best exposure time and film to focus distance (FED) 
for each specimen thickness. The X-ray source used was 
capable of varying both the voltage and the current which 
produce the X-rays. It was quickly discovered that each of 
the variable parameters must be chosen correctly if an 
acceptable degree of resolution was to be realized. Many 
test shots were taken of a given laminate at various cur- 
rents, voltages, distances, and exposure times. The nega- 
tives from these tests were then used to make radiographic 
prints to determine the correct exposure parameters. 

It was found that the single-ply and the two-ply speci- 
mens could be x-rayed at an FED of 50 cm with a current of 5 
mA. On the thicker four-ply and eight-ply specimens, the 
best radiographs were taken at an FED of 80 cm with a cur- 
rent of 5 mA. The higher FFD serves to place the X-ray 
focal point further away from the object, thus making the 
X-ray seem more like a point source; that is, the rays are 
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almost parallel when they strike the specimen. The voltages 
and exposure times used varied depending on the thickness of 
the panel. A graph of current-exposure time product versus 
the specimen thickness for various voltages is shown in Fig- 
ure (2). Any X-ray exposure taken with data from this graph 
was found to give excellent resolution on the radiographic 
print. As an example, the four-ply panels, with a thickness 
of 0.712 mm, were x-rayed for three minutes at 30 kV with a 
tube current of 5 mA, giving the required product of fifteen 
mA-min as determined from the graph. 

The image produced from the X-ray is captured by simply 
securing a sheet of Type 55 Polaroid film to the back of the 
specimen at various load levels during the test. It must be 
emphasized that the film is not placed in a film holder, 
such as a Polaroid Model 455, which has an intensifying 
screen. The intensifying screen is not needed for this type 
of exposure; in fact, use of the film holder actually made 
the resolution worse on the radiograph. Thus the sheet of 
film is taped to the specimen encased only in the light 
paper provided by the manufacturer. 

The number of X-rays absorbed or transmitted by a mate- 
rial is dependent on its atomic number. For example, an 
X-ray of a graphite/epoxy panel will produce no distinct 
fiber pattern on a radiograph since both the fibers and the 
matrix are composed of essentially carbon. In fact, boron 
(atomic number = 5) and aluminum (atomic number = 13) 


CURRENT-EXPOSURE TIME PRODUCT (niA-min) 


OF 


^ it .i* 


QUALVl 



Figure 2. X-ray Exposure Chart. 
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transmit amounts of X-rays that are so close that it is dif- 
ficult to distinguish between the two of them on a radio- 
graph. The boron fibers commonly available, however, are 
manufactured by depositing the boron on a thin tungsten 
filament, and this tungsten (atomic number = 74) core will 
appear as a very dark black line on a radiograph. Thus the 
tungsten core of each fiber can be easily identified on the 
radiographic print. 

A typical radiograph of a single-ply boron/aluminum 
panel is shown in Figure (3). The dark lines are the tung- 
sten cores of the boron fibers, while the light area around 
these cores is the boron itself. The medium gray area 
between the fibers is the aluminum matrix. Thus as the 
specimen is loaded and fibers begin to break, it is a simple 
matter to count the number of broken fibers ahead of the 
notch at various load levels. It is also possible to obtain 
a very accurate measure of the crack opening displacement 
(COD) by merely scaling the appropriate fiber displacements 
directly from the radiographic print made in the darkroom. 
COD measurements were also taken with a clip gage which is 
specifically designed to measure crack opening displacement 
for notched specimens. The output of this gage can be care- 
fully calibrated such that the strain output from the gage 
can be converted to a displacement in millimeters. 
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Radiograph of single-ply boron/aluniram composite panel. 
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Brittle Coating Technique 

A brittle lacquer was applied to each specimen to meas- 
ure the amount of longitudinal yielding between the last 
broken fiber and the first unbroken fiber at the end of the 
notch. The brittle lacquer used in these tests was the com- 
mercially available Tenslac, manufactured by the Measure- 
ments Group in Raleigh, North Carolina. This lacquer will 
crack when a certain threshold strain is reached in the 
boron/aluminum panel. This threshold strain can be affected 
by the coating thickness as well as the temperature and 
humidity of the test environment, 

As previously mentioned, each specimen was carefully 
cleaned and polished to provide a specimen surface that was 
free from surface flaws. The polished surface also provided 
a highly reflective surface below the brittle lacquer. This 
made the cracks reflect the light much more clearly so that 
the extent of cracks could be observed and photographed. It 
was noted that the ability to see a given crack in the coat- 
ing was strongly affected by the angle at which the light 
approached the surface of the specimen. The cracks could 
not be seen if the light approached the surface along a line 
normal to that surface. Rather, the best results were 
obtained when the approaching light made an angle of approx- 
imately 30 degrees with the surface normal. 

The brittle lacquer was applied to each specimen in up 
to twenty very thin coats to minimize bubbling and provide 
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the desired "orange-peel" finish. The final thickness of 
the lacquer coating was approximately 0.075 millimeters 
which gave a threshold strain of 1200 to 1800 microstrain 
depending on test conditions. The threshold strain for each 
specimen tested was measured by using strain gages that were 
mounted on each panel on the opposite side from the coating. 
The remote strain at which the cracks appeared directly 
opposite the gages was taken as the threshold strain for 
that specimen. This procedure eliminated the need of pre- 
paring calibration bars or even assuming a nominal threshold 
strain. Rather, the exact threshold strain at the time of 
the test for each specimen could be accurately measured and 
recorded. The brittle lacquer coating also gave a very 
clear picture of the stress state at the remote sections of 
the specimens. If the remote stress was truly uniform, the 
cracks in the coating would be straight and perpendicular to 
the fibers outside the region affected by the defect. The 
output from the strain gages was measured through the use of 
a Vishay VE-20 Strain Indicator and a Vishay VE-21 Switching 
and Balancing Unit. 
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General Testing Procedure 

After allowing the brittle coating to dry for 18 to 24 
hours, the specimens were loaded statically with a Baldwin 
tensile testing machine. A photograph of the test set-up is 
shown in Figure (4). A baseline X-ray was taken of each 
panel before the loading cycle was begun so that any damage 
caused by machining or polishing could be noted. Then as 
the loading sequence began, photographs were taken of the 
cracks in the coating at various load intervals with a Nikon 
FM camera equipped with a close-up lens. A bright light was 
used to illuminate the crack pattern so that a small lens 
aperture could be set to increase depth-of-field. 

X-rays were taken in the manner described earlier at 
selected load levels so that the number of broken fibers in 
the specimen could be determined. The majority of the X-ray 
exposures were taken in the last half of the load cycle 
since the first fiber breaks occurred at approximately 50 
percent of the failure load. The crack opening displacement 
was measured both from the radiographs and with a clip gage 
specifically designed to measure COD. 

A total of ten to twelve X-rays and twelve to fifteen 
brittle coating photographs were made during each test. The 
strain output from each gage was recorded at numerous load 
levels as were the load readings from the Baldwin machine. 
This testing procedure was repeated for boron/aluminum pan- 
els of many different widths, thicknesses, and initial notch 


widths . 
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Figure 4. Experimental Test Set-up. 


CHAPTER III 


ANALYSIS 

Mathematical Model Description 
Though the main emphasis of this paper is the experi- 
mental investigation, it is important that the mathematical 
model which motivated the present study and the assumptions 
made therein are clearly understood. Thus at this point, a 
brief description of the equations involved in the deriva- 
tion of the shear-lag model will be presented. Following 
Goree and Gross [1], the laminate is modeled as a two-dimen- 
sional region, having a single row of parallel, identical, 
equally spaced fibers, separated by matrix as shown in Fig- 
ure (5). The initial damage is taken to consist of an arbi- 
trary number of broken fibers such that all breaks lie along 
the x-axis. The longitudinal matrix damage is introduced at 
the end of the initial notch as in [1]. The additional 
notch tip transverse damage consists of an arbitrary number 
of broken fibers which are constrained by the adjoining 
matrix and/or the unbroken fibers through the thickness. 
These fibers in the transverse damage zone will be referred 
to as constrained fibers. It is mathematically untractable 
to account precisely for the distribution of fiber damage 
occurring in the region ahead of the notch. The model 
assumes that all the breaks occur on the x-axis and accounts 
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Figure 5. Two -dimensional unidirectional lamina with broken 

fibers, longitudinal matrix splitting and yielding, 
and transverse fiber and matrix damage (first 
quadrant) . 









for the stiffness by assuming that these fibers carry a 
reduced load by defining a stiffness coefficient, X, such 
that 
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y - stress in the constrained fiber 
“ stress in the first unbroken fiber ' 

The free-body diagram shown in Figure (6) is considered 
with additional conditions for the last broken fiber, 
denoted by n = N, and y < L, to account for the longitudinal 
damage taken as 


x 


*N+1 


t <y - £>, 

o J 


( 1 ) 


where 


<y - z> “ 1, y k Z, and 

<y - z> = 0, y < z. (2) 

L is defined as the total damaged length, Z as the split 
length, and t as the matrix yield stress. Yielding is 
assumed to occur when the matrix reaches the yield strain, 

. Splitting occurs at a multiple of X Q as given by the 
particular matrix material. 

The equilibrium equations in the longitudinal and 
transverse directions for all fibers n, with the exception 
of N and N + 1 when y < L, are then 
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Figure 6. Free-body diagram of a typical element. 
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( 3 ) 


(4) 


For fiber N, when y i L, t In+1 s “ t 0 <y “ l> ' and tlie 
equilibrium equations are 


Ap ^°fIn 
T cly 


V y 


fc> 



(5) 


and 


a I - a L + ^ <y - i> + x| N > = 0. (6) 

m'n+1 m'N 2Hy o 'N 


For fiber N + 1/ when y £ L> t| n+ -^= 
the equilibrium equations are 


t q <y - <t>, and 


a F + ,1 + t <y “ «.» = 0 , 

t ay N+l o 


(7) 


and 


°M ! N+2 ' 0 M* N+l + ?"3y {t N+2 " T o 


h d <tU., - <„<y - «» = °- (8) 


The stress-displacement relations common to the shear- 


lag assumption, 
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(9) 

( 10 ) 
( 11 ) 


are introduced into the equilibrium equations, leading to 
the replacement o£ the three stress quantities by the axial 
and transverse displacements. These assumptions, however, 
simplify the equilibrium equations by removing the tran- 
sverse displacement dependence from the longitudinal equi- 
librium equation. The fiber stress and the matrix shear 
stress can then be determined without solving the trans- 
verse equilibrium equation. Only the longitudinal equilib- 
rium equation will now be considered. For all fibers, 
except N and N+l where y S L, the longitudinal equilibrium 
equation becomes 


E F A F h d v n 


G M t 


dy' 


+ v . . - 2v + v n = 0, 
n+l n n-1 


( 12 ) 


For fiber N where y £ L, 

_ v 

V dy 2 «- 1 N 



t <y — (t> 
o 


0. 


(13) 


For fiber N+l where y £ L, 


EpAph d 2 v N+1 ^ ^ 

"ST"d?“ N+2 


V N+1 



t <y 
o 


i> 


0. (14) 
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It is convenient to use the non-dimensional forms of 
Equations (12), (13), and (14). Noting the coefficient of 
the second derivative term in the above equations, define 


y = 


V 


n and o 


F'n 


= oo 


n 


The other required quantities, such as the axial displace- 
ment (v ) and the matrix yield stress (i ), can also be 

n o 

written in their non-dimensional forms. From the defini- 
tions o* y and o r I it follows that 
■* r ' n 
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The quantities n, t, o n , V n , a, and & are non-dimensional, 
whereas Ep , Ap , t, L, and t represent the actual fiber modu- 
lus, fiber area, lamina thickness, and damage zone dimen- 
sions, respectively. 


The resulting equilibrium equations in non-dimensional 
form are given by 


d 2 V 


dn 


r* V n+1 ' 2V n * V n-1 = °' 


d 2 V, 


dn 

d 2 V 


r * v n + i - 2V n * Vi * - f< ">' * nd 


dn 


* V N+2 - 2 Vl * V N * £ <"> 


(16) 


(17) 


(18) 


where f(n) is a new unknown function such that 


f(n) = V N - V N+1 - t q <n - &>, n < o, and 

f(n) =0, n * a. (19) 


These differential-difference equations may be reduced to 
differential equations by introducing the even valued trans- 
form as 


V(n,0) = V n (M/ 2 ♦ I V (n) cos(n0 ) , (20) 

0 n=l n 

from which 


2 — 

V (n) = - J V(n,0) cos ( n0 )d0 . (21) 

n v Q 

Making use of the above transformation and the orthogo- 
nality property of the circular functions, the three 
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or 
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.1 
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equilibrium equations may be written as one equation valid 
for all values of n and n as 


2 " j2y 2 

Z / { -211 - cos(9))V) cos(nO )d 0 ■ - X 

71 0 dn 


u 

<o-n> / f (n)(cos( (N^l ) e 1 

0 


cos(NO)) cos(n0)d0. 

( 22 ) 


This equation is of the form 

2 v 

— J F(n, 6 ) cos(n 0 )d 0 * 0 for all n and n. 

’ 0 

Noting the definition of V(n,0) in Equations (20) and 
(21) it is seen that the function F(n,0) is even valued in 0 
and therefore, if the integrand is to vanish for all n, the 
function F(»t,0) must be zero. The single equation specify- 
ing V(n,0) is then 

2 — 

^-7 - 6 2 V = - <o-n> D 2 f ( n ) , (23) 

dr/ 


2 2 

6 = 2(1 — cos( 0 )) = 4 sin (0/2), and 

2 

D = cos(NO) - cos[(N+l)0]. 


where 
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The solution to the problem of vanishing stresses and 
displacements at infinity and uniform compression on the 
ends of the broken fibers will now be sought. The complete 
solution is obtained by adding the results corresponding to 
uniform axial stress and no broken fibers to the following 
solution. The appropriate boundary conditions are thus 


V n = 0 


(24) 


Tr T = 0 as 


and 


(25) 


V n = 0 for n = 0, for unbroken fibers. 


(26) 


n — 
SW = °n = 


= -1, for broken fibers (n=0 to n=N), (27) 


n M+l 

— = -1 ♦ j — 3ri ' ^ or constrained fibers 

(n=N+l to n=M). 


(28) 


The complete solution to Equation (23) satisfying van- 
ishing stresses and displacements at infinity is given by 


V( n , 0 ) = A(0) e 


J (sinh(6(r;-t) )<o-t> f(t)dt. 


(29) 


where the unknown functions are A(0) and f(t). The remain- 
ing boundary conditions give 
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dV n (0) 

dn 


« IT n Ot 

z I { -6A( 0 ) ♦ D / cosh(0t)f (t)dt)cos(n8)d0 = -1, 
11 0 0 

n = 1, 2, . . . , N, (30) 


dV (0) o v 2 a 

— 4- — =£ J { -6A( 0 ) ♦ D / cosh(6t)f (t)dt) X 
an ti o 0 


cos(n0)d0 = -1 ♦ ,0 f*M+ 1' n “ N+ *' •••* 

(31) 


2 71 

V n (0) = ± / { A( 0 ) 
n v 0 



a 

J sinh(6t)f (t)dt) cos(n0)d0 = 0, 

0 

n = M+l, (32) 


Equation (32) is solved exactly by taking 


2 a 00 

A(0) — t- I sinh(6t)f ( t)dt = I B_ cos(m0) 

6 0 m=0 m 


(33) 


where m is a broken fiber index and the B_ are constants. 

m 

The stress in the first unbroken fiber, °f!m+ 1 g * ven by 


d Vl 2 " 


a 


-n— = - =- J {-6A( 0 ) + D / cosh( fit ) f ( t )dt) X 

an IT 0 0 


cos( (M+l )0 ]d0 = -1 ♦ °f!m+1' 


(34) 


Using Equatio: s (33) and (34) in Equations (30) and 
(31), A(0) and may be eliminated, resulting in M+l 
algebraic equations for M+l constants B m in terms of f(n) 


which is, as yet, unknown. For longitudinal matrix damage. 
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Equations (30) and (31) must be supplemented by the condi- 
tions that 

f ( n ) = g(n) - T <n - &>, n < a, and (35) 

g(o) = T , (36) 

where 


9(i) 



V 


N+l ' 


From Equations (28) and (33), A(6) may be eliminated to 
obtain V(n,0) in terms of constants B m and the unknown func- 
tion f(t). Recalling the relation between V(ti, 8) and V n (n), 
an expression can be obtained for the axial fiber displace- 
ment v n ( T >) as 


V> 


0 V <r M 

= ± / e"° n r B cos(m0)cos(n0)d0 
v 0 m=0 m 



a 

/ f(t) {c ( | t-n I ) - c ( t+ri ) ) dt, 

.. n n 


where 

c n (S) = I y 1C e ~ ^ cos(n0)d0. 


( 37 ) 
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Equations (30) and (31) then become 
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- t / e dt) cos(n0)d0 = -1, 
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n = 0, 1, . . . , N, 


(38) 


and 

1 ♦ l, 
(39) 

Equations (36) and (37) give 
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which is a Fredholm integral equation of the second kind. 

The last condition that must be satisfied is that 

g(o) = T. (41) 

The series equation and the integral equation are cou- 
pled by the Fourier coefficients and the function g(n) and 
may be solved on a digital computer using a Gauss-Legendre 
quadrature scheme. Recent improvements in the solution 
scheme have reduced the computation time to less than two 
minutes, even for a very large number of broken fibers and 
large matrix damage. 

Brittle Coating Interpretation 

The mathematical model described in the previous sec- 
tion gives the fiber displacements and stresses throughout 
the laminate as well as the longitudinal yield region 
between the last broken fiber and the first unbroken fiber. 
An important step in relating the behavior predicted by the 
model to the measured laboratory values is thus to compare 
the measured longitudinal yield zone with the predicted 
value of the yield length from the model. The difficulty in 
this comparison is the interpretation of the crack pattern 
found in tne brittle coating photographs. 

Using the fiber stresses and displacements from the 
mathematical model it is possible to approximate the surface 
stress state in the aluminum matrix between the fibers. 
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Since the thin coating of brittle lacquer is sensitive to 
strains in this small region, the angle of the cracks in the 
coating should correspond to the principal strain (or 
stress) angles calculated from a simple Mohr's circle analy- 
sis for the aluminum. It should be noted that the matrix 
normal stress in the vertical direction was omitted in the 
shear-lag analysis for the fiber stresses since this stress 
is small compared to the fiber stress. However, it is 
appropriate to now consider this matrix stress since it is 
part of the stress state in the matrix which influences the 
crack pattern in the coating. 

A good approximation for the average matrix normal 
strain in the vertical direction can be calculated by aver- 
aging the strains in the two adjacent fibers. 


, dV n + dV n+l , 
e y n ’ 21^ ( “3TT <Jn ] ' 


(42) 


The horizontal matrix strain is assumed due to a Pois- 
son contraction alone, and is thus derived by multiplying 
Equation (42) by the Poisson ratio for aluminum, 0.3. The 
stresses in the matrix are then and o y = 0. 

Finally, the matrix shear stress is calculated from Equation 
( 15 ) as 


T = (V 

n — n 


w 


T 


(43) 
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Thus a typical element as shown in Figure (7) r>elow has 
a known state of stress and the angle through which the ele- 
ment must rotate to reach the angle of principal stress can 
be calculated from 


tan( 26 ) 


(44) 


Figure 7. Mohr s Circle Analysis. 


The growth of longitudinal damage predicted from the 
above analysis can now be compared to the growth from actual 
specimen photographs. The longitudinal crack pattern for a 
specimen having an overall width of 73.2 mm and an initial 
notch width of 12.7 mm will be examined at various load lev- 
els and compared to predicted strain data from the model. 

The first load step examined is very early in the cycle 
before any longitudinal yielding has begun. At an axial 
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load of 800 N the model predicted that very few cracks 
should form, since the strain falls below the threshold 
strain of 1500 we at a y-value of approximately 1 mm. The 
model further predicts that these first cracks that form in 
the coating will form at an angle of 37 degrees from the 
horizontal. Table (II) shows the predicted data from the 
model at the load levels to be examined in this analysis. 
Figure (8) shows the appearance of the coating at a load of 
800 N. It is clear that these first cracks, forming an 
angle measured to be 32 degrees, appear directly above the 
end of the notch between the last broken fiber and the first 
unbroken fiber. Moreover, the uppermost crack forms at a 
y-value on the order of that predicted by the strain values 
computed from the analysis. 

At a load of 3000 N the model predicts that approxi- 
mately 2.5 mm of longitudinal matrix will yield and that the 
principal strain in the matrix will exceed the brittle lac- 
quer threshold strain value for y-values less than 5 mm, 
with the last predicted crack forming at an angle of 30 
degrees. Examination of Figure (9) reveals that the upper- 
most crack, at a y-value of 6 mm, forms an angle of 29 
degrees with the horizontal. In this figure it can also be 
seen that a region of very intense cracking has formed, con- 
sisting of short, closely spaced cracks and lacquer crazing. 
The length of the region containing these short cracks cor- 
responds very closely to the yield length of 2.5 mm pre- 
dicted by the model. 
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The length of the longitudinal matrix yield region at 
4S00 N is predicted by the model to be approximately 5 mm, 
and the uppermost crack in the coating at this load should 
form 8 mm above the notch at an angle of 29 degrees. Figure 
(10) shows that the uppermost crack is present 12 mm above 
the notch at an angle of 27 degrees and that the region of 
severe cracking and crazing has grown to approximately 5 mm. 

Finally, in Figure (11) the uppermost crack forms 14 mm 
above the notch at an angle of 20 degrees with the horizon- 
tal, and the region of severe cracking has grown to 10 mm at 
this load level of 6500 N. The predicted location and 
orientation of the last crack at this load level are 15 mm 
above the notch at an angle of 22 degrees with the horizon- 
tal. The yield length predicted by the mathematical model 
at this load level is 10 mm. Thus the stresses and strains 
in the matrix actually measured in the laboratory appear to 
agree reasonably well with those predicted by the mathemati- 
cal model. The growth of the crack region as well as the 
severely cracked region is shown in graphical form in Figure 
(12). It is clear that the brittle coating provides a very 
good qualitative picture as well as some quantitative data 
concerning the yield region at the end of the notch and its 
growth as the load is increased. 

It should be noted that the cracking discussed above 
occurs relatively early in the load cycle before the first 
fiber break has occurred. As fibers begin to break in the 


43 


laminate, the loading is stopped and restarted each time an 
X-ray exposure is made. Ar is the case with all brittle 
lacquers, the load must be r« moved completely each time the 
loading is stopped if the threshold is to remain valid. Of 
course, this would greatly alter the stiffness of the previ- 
ously yielded aluminum matrix and thus change the behavior 
of the laminate. The loading cycle was thus continued after 
each X-ray exposure without first removing the load. 

The brittle coating, while providing actual numerical 
approximations in the early part of the load cycle before 
the first fiber has broken, can provide only qualitative 
results later in the cycle. It is felt, however, that the 
severely cracked region is closely associated with yielding 
in the matrix. Several tests were performed usit.g pure alu- 
minum notched panels and in each case the region(s) in which 
yielding was known to occur exhibited this same intense 
cracking. The region of intense cracking and flaking is 
thus used as indicator of matrix yielding and will be refer- 
red to as the severe damage length in the remainder of this 


paper . 
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Figure 8. Brittle coating crack pattern for an axial load 
of 800 N. 
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Brittle coating interpretation data in graphical form 


CHAPTER IV 


RESULTS 

Damage Growth Sequence 

The main objective of this paper is to determine accu- 
rately the properties of boron/aluminum laminates and then 
compare the predicted values to those measured in the labo- 
ratory using actual specimens. The experimental tests 
should provide values for the matrix yield length, all fiber 
displacements, fiber stresses, and stable notch extension. 
Material properties and laminate dimensions from handbooks 
and manufacturer's specifications should allow the experi- 
menters to convert the non-dimensional predicted values from 
the model to real values suitable for comparison with expe- 
rimental results. 

As previously discussed, the use of radiographs and 
brittle coating photographs will present a very detailed 
picture of the growth of damage during a loading cycle of 
any unidirectional boron/aluminum laminate. It is therefore 
appropriate to begin the examination of the results of this 
study with a step-by-step presentation of a typical set of 
photographs of one of the experimental specimens. The spec- 
imen chosen for this analysis is typical of all of the test 
panels, exhibiting both matrix yielding and stable notch 


extension in the form of broken fibers ahead of the notch. 
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The set of photographs and radiographs show the pro- 
gression of damage in a two-ply boron/aluminum laminate 
which is 73.1 mm wide. A 19.05 mm wide notch (109 broken 
fibers) was cut perpendicular to the fiber direction in the 
center of the panel before the testing procedure was begun. 
Brittle coating photographs and X-ray exposures were made at 
selected intervals in the loading sequence. The series of 
radiographs and photographs is presented in Figures (13) 
through (28) beginning at the end of this section. 

As previously mentioned, a baseline X-ray was taken of 
each specimen before the loading cycle was begun to detect 
any damage which may have been caused by machining or pol- 
ishing. The baseline X-ray for the two-ply specimen is 
shown in Figure (13). The right end of the notch can be 
seen in this radiograph, along with the fiber cores which 
run perpendicular to the notch. The initial notch opening 
was 0.381 mm; thus crack opening displacement at each load 
level could be measured from the radiographs once the scale 
factor had been determined. It is clear that there are more 
fibers visible in this radiograph than in the earlier sin- 
gle-ply radiograph. This indicates that the fibers in each 
ply, as expected, are not perfectly aligned with the fibers 
in adjacent plies. It can be seen, however, that there are 
usually two fiber cores very close together on the radio- 
graph, indicating the two fibers in adjacent plies. 
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In the first few photographs, the cracks in the brittle 
coating grow in the longitudinal direction, even though no 
fiber breaks have occurred. In Figure (17) it can be seen 
that the first fiber breaks have occurred, indicated by the 
breaks in the fiber cores near the end of the notch. It 
should be noted that from this point there will be no longi- 
tudinal growth of the severe damage region, since the first 
previously unbroken fiber is now free to displace. The pro- 
gression of damage, indicated by severe cracking and flaking 
of the brittle lacquer, will now be in the transverse direc- 
tion as more fibers continue to break. 

In Figure (24) the severe damage region has grown tran- 
sversely from the notch, clearly indicating the existence of 
more than one longitudinal damage region in the matrix. It 
is this phenomenon that allows the COD to grow at a faster 
rate than predicted by the model, since the model at the 
present time can include only one yield region. The last 
figures in the damage sequence indicate that the transverse 
damage region continues to grow until the failure load of 
the laminate is reached. It should be noted that the exten- 
sion of the notch in the form of broken fibers does not 
occur in a straight line parallel to the notch. Rather, the 
crack grows in a jagged fashion with tbs fiber breaks occur- 
ring at different vertical locations as the crack grows and 
more fiber breaks are observed. The fibers in this region 
are thus broken but still able to carry some reduced amount 





original PMW h 

BUCK AND VW.TE PHOTOGRAPH 


54 







Figure 13. Baseline radiograph of a notched, two-ply boron /aluminum composite panel. 
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Figure 14. 


Brittle lacquer crack pattern at 2670 N showing 
formation of first cracks in the coating. 
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Figure 20. Brittle lacquer crack pattern at 9350 N showing 
transverse growth of severe damage region. 
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Figure 26. Brittle lacquer crack pattern at 12,800 N. 
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Brittle lacquer crack pattern at 13,200 N 
showing large transverse progression of the 
severe damage region. 
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0»nT>l Results a nd Discussion 

By examining the sequence of photographs and radio- 
graphs previously discussed, one can form a very clear pic- 
ture of the manner in which damage accumulates in a unidi- 
rectional laminate as the failure load is approached. 

Similar sequences could be presented for the entire range of 
specimens since the same testing procedure was followed for 
every specimen used in the study. The two-ply sequence 
presented, however, was quite typical of all of the tests 
and it was thus felt that the remainder of the results could 
be presented more effectively in the form of graphs which 
show results of all of the tests. Each of the questions 
raised in the introduction will be answered and discussed 
thoroughly in the remainder of this section. 

The question of thickness effect on the behavior of a 
unidirectional laminate can now be answered by examining the 
results of each set of tests. A full range of specimen 
thicknesses was tested for each panel width and initial 
notch width. It was observed that the severe damage region 
as well as the overall appearance of the brittle coating at 
a given remote strain was very similar for the single-ply, 
the two-ply, the four-ply, and the eight-ply. The length of 
the severe damage region versus the number of plies is plot- 
ted for various numbers of crack lengths in Figure (29). 
Excellent agreement between the length of the observed dam- 
age region and that predicted by the model exists as dis- 
cussed in the previous section. 


71 


It was also noted that the dimensions of the transverse 
damage region were very similar for a given initial notch 
width. The number of broken fibers ahead of the notch was 
approximately the same for the one-, two-, four-, and 
eight-ply laminates, as one can see in Figure (30). Thus, 
from these two graphs, is it clear that the longitudinal 
yielding and transverse damage regions depend very little on 
the number of plys in the unidirectional laminate. Again, 
the number of fibers broken ahead of the notch agrees well 
with the number of constrained fibers predicted by the math- 
ematical model. 

Figure (31) shows the crack opening displacement (COD) 
for the set of specimens with an initial notch of nineteen 
broken fibers (or 3.17 mm). It is clear that the COD for 
each thickness is almost the same for a given percentage of 
the failure load and therefore the same remote strain. It 
can also be noted that the first fiber break occurs between 
50 and 55 percent of the ultimate failure load for each 
thickness tested. This phenomenon was observed in all of 
the tests regardless of notch width or specimen thickness. 
Figure (31) also indicates that the maximum COD at failure 
was almost identical for the four thicknesses tested. 

Figures (32) and (33) show the COD behavior for speci- 
mens with 53 and 89 initially broken fibers, respectively. 

It is again clear that the specimen thickness has little 
effect on the COD for a given percentage of the failure 
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load. The first fiber breaks occur at roughly the same 
remote strain level for a given notch width, as shovm in the 
earlier graph. The appearance of the brittle lacquer coat- 
ing was observed to be the same for each thickness as in 
earlier tests. Also, the stress necessary to break the 
first fiber, the length of the severe damage region, and the 
notched strength of each laminate all are in close agreement 
with the quantities predicted by the mathematical model. 

This agreement will be discussed in detail later in this 
section. The similarity in behavior regardless of the spec- 
imen thickness (for up to eight plies) was observed for all 
initial notch widths, thus clearly indicating that a thick 
panel with unidirectional fibers in each ply can indeed be 
modeled as an equivr lent monolayer. 

Another important objective of the present study was to 
determine the effect of different matrix properties on the 
fracture behavior of unidirectional laminates. Toward this 
end, certain sets of panels (as detailed in Table I) were 
heat-treated to bring their matrix properties close to the 
6061-T6 condition instead of the fully annealed 6061-0 con- 
dition as in the other panels. It was felt that a matrix 
with a different effective stiffness and yield strength 
would alter the fracture behavior since it is the matrix 
alone that transfers load to adjacent fibers in the lami- 
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Figure (34) clearly shows that again the COD at a given 
percentage of the ultimate failure load is not strongly 
dependent on the number of plys in the laminate for the 
6061-T6 panels. As in the 6061-0 tests, the first fiber 
break occurs at SC o 55 percent of the failure load. The 
COD at failure for -he four thicknesses is almost the same, 
though it should be noted that this COD is much less than 
the COD for a similar panel with an annealed aluminum 
matrix. This change will be discussed fully, as it will 
help to determine the effect of matrix properties on the 
damage growth. 

Figure (35) is a remote stress-COD curve for various 
notch widths in a four-ply, 6061-0 matrix panel. As 
expected, specimens with short notches are stiffer than 

i 

specimens with longer notches. Figure (36) shows the COD 
behavior for the same nr tch widths in a four-ply, 6061-T6 
Matrix panel. It is clear that there is a drastic differ- 
ence in the notch opening under load when the properties of 
the matrix are changed. The stiffer and stronger 6061-T6 
matrix does not allow the fibers to displace nearly as much 
as the fully annealed 6061-0 matrix. 

It was suggested in (2) that the mathematical model 
employing the shear-lag assumption could predict accurately 
the load-COD curves for specimens of various notch widths if 
certain laminate material properties were known. Specifi- 
cally, these are the normalized matrix yield stress and the 
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effective shear stiffness of the fiber-matrix region which 
are calculated from an experimental load-COD curve. By 
using these properties, the load-COD behavior for any number 
of broken fibers could then be predicted accurately. Agree- 
ment was demonstrated in (2) for the narrow specimens and 
short notch lengths presented by Awerbuch and Hahn |6). 

Using the much more extensive results of the present study 
it can be shown that, until large transverse damage is pres- 
ent, the load-COD curves for linger notch widths can also be 
matched. 

In Figure (37) a typical load-COD curve for a 4-ply 
unidirectional boron/aluminum composite panel is given for 
specimens having an overall width of 73.1 mm and 69 ini- 
tially broken fibers. Both a specimen having a S061-0 ann- 
ealed matrix and a 6061-T6 matrix are shown, along with the 
load-COD curves predicted by the model. Using the values 
for the yield stress, and the effective shear stiffness, 

G^/h, given in |2|, the load-COD curve for the 6061-0 speci- 
men can be closely matched at lower load levels even for 
longer notches than in |6|. The agreement is not as good 
for large notch extension as discussed later in this sec- 
tion. It must now be determined if appropriate changes in 
the two properties mentioned above will allow the model to 
predict accurately the load-COD curve for the 6061-T6 speci- 
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The load-CGD curve t*>r the 6061-T6 panel can be matched 

by increasing the value of the matrix yield stress 3.5 times 

2 

to 308 X 10 N/m and the value of the effective shear stiff- 

3 

ness 1.5 times to 172.5 X 10 N/m . These changes not only 
seem qu? te reasonable but the yield stress change agrees 
well with standard aluminum properties. The change in the 
value for G^/h is attributed more to a change in the shear 
transfer length, h, than to a change in which should 
remain relatively constant with any heat-treatment . The 
yield stress for 6061-T6 aluminum, on the other hand, is 
roughly 3-4 times the yield stress of 6061-0 aluminum. 

The load-COD curve predicted for the 6061-T6 specimen, 
in fact, matches the actual one better than the curve pre- 
dicted for the 6061-0 specimen. It must be noted that both 
curves fit quite well at lower load levels; at greater load 
levels, however, the model predicts a higher laminate stiff- 
ness than is observed on the actual specimen in the labora- 
tory. The deviation from the predicted curves occurs after 
several fibers have broken and the notch has begun to extend 
a significant amount. Since it was observed in all cases 
that the notches in the specimens having annealed matrices 
ran much further before fracture than their 6061-T6 counter- 
parts, the discrepancy can easily be explained. It is the 
matrix yielding and broken fibers that allow the notch to 
open, thus the COD for the annealed matrix specimens ’..’ill be 
considerably greater than tho COD for the T6 matrix at a 
load level corresponding to significant damage. 
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It was mentioned in the previous section that the 
appearance of the brittle coating indicated that there may 
be more than one yield region near the end of the notch. As 
each successive fiber break occurs, the newly broken fiber 
can move with respect to the unbroken adjacent fiber. Thus 
a new region of significant displacement and matrix yielding 
will form, allowing the COD to increase more than if the 
matrix remains elastic. At the present time the models of 
\ 2 ) can incorporate only one yield region between the last 
broken fiber and the first unbroken fiber. Thus, if more 
than one yield region forms, the model cannot predict accu- 
rately the displacements of fibers after the notch extension 
has begun. 

Future plans include the modification of the mathemati- 
cal model to include a new yield region as each successive 
fiber is broken during a loading cycle. However, substan- 
tial information can be obtained by summarizing the present 
experimental data. Figure (38) shows the crack opening dis- 
placement at failure as a function of the number of ini- 
tially broken fibers. Clearly the crack opening displace- 
ment for the 6061-T6 panels is substantially less than that 
of the fully annealed 6061-0 panels. In fact, the differ- 
ence becomes larger as the number of initially broken fibers 
is increased. This phenomenon is caused by the relatively 
large notch extension of tne fully annealed specimens. The 
stronger T6 matrix cannot absorb as much energy through 
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yielding as the annealed matrix, thus the load is transfer- 
red to the next unbroken fiber sooner and the notch exten- 
sion becomes unstable, leading to earlier failure. 

The longitudinal damage length as described in an ear- 
lier section is plotted as a function of the number of bro- 
ken fibers in Figure (39). The predicted value for this 
damage length is obtained by calculating the yield length 
necessary to reach the breaking stress in the next unbroken 
fiber through the use of the mathematical model. This value 
should correspond to the severe damage region on the experi- 
mental specimens since there should be no further yielding 
in the region between fibers once the next fiber has broken 
and is free to displace. It is clear that the predicted 
trend of damage growth from the model is a very good indica- 
tor of the actual damage present in a unidirectional lami- 
nate. The stronger T6 matrix substantially reduces the 
yielding in the matrix and thus reduces the amount of damage 
tolerance possessed by a boron/aluminum laminate. The 
result of increasing the yield strength of the matrix mate- 
rial used in a unidirectional laminate is thus to make the 
laminate more notch sensitive and actually weaker than a 
laminate having a more ductile, energy-absorbing matrix. 

The notched strengths of the test panels and the reduc- 
tion in strength for the 6061-T6 matrix can be studied in 
Figure (40). It should again be noted that the effect of 
the change in matrix properties becomes much more evident as 
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the notch becomes larger; that is, the number of initially 
broken fibers is increased. As in (2], the strength curve 
for the fully-annealed 6061-0 matrix specimens can be pre- 
dicted by the mathematical model using 7 constrained fibers 
ahead of the notch. It should be remembered that in the 
damage photograph sequence shown earlier there were approxi- 
mately 7 constrained fibers ahead of the notch tip even for 
large notch extension. The strength curve for the 6061-T6 
panels can be approximated by using 3 constrained fibers 
ahead of the notch. This change has physical significance; 
the stable notch extensicn is larger for the fully annealed 
specimens than for the heat-treated specimens. Thus as dis- 
cussed earlier, the heat-treating of the matrix actually 
reduces the damage tolerance capabilities of the laminate 
and makes the unidirectional laminate weaker than one with a 
fully annealed matrix. 

One should be careful in extending this reasoning to 
even stronger matrix materials in which the transition is 
made from longitudinal yielding to longitudinal splitting. 
For example, an epoxy matrix is one in which little yielding 
is exhibited but longitudinal splits occur at the notch tips 
under low load. This phenomenon is discussed in (2) and has 
been observed by many investigators. It is important to 
realize that once this small split, forms, it will quickly 
grow to the machine grips in an experimental test (or to 
infinity as predicted by the model). The stress 


concentrations are then removed and the notched laminate 
will fail at the net section failure load. Thus the conclu- 
sion that a stiffer or stronger matrix will weaken the 
notched laminate is valid only for mf.trices that are not 
brittle enough to form longitudinal splits between fibers at 
the end of the notch. Some small splitting was observed in 
the present study in 6061-T6 laminates with very long 
notches, possibly indicating that the T6 matrix is close to 
the transition from matrix yielding to splitting. 
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Figure 29. Variation of the severe damage length with the 
number of plies for selected numbers of broken 
f ibers . 
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Figure 30. Variation of the transverse notch extension 

with the number of plies for selected numbers 
of broken fibers. 







100 


ORfGWAt PAGE IS 
OF POOR QUALITY 



3 


to 

I— * 

Q 

O 


x: 

UJ 

a, 

o 

XL 

CJ 

2 

CJ 


avo'i 3am i vi imoaad 


Figure 33. Percent failure load versus crack opening displacement for 89 broken fibers 
per ply; comparison of various numbers of plies for a 6061-0 Al matrix. 
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Figure 36. Remote stress versus crack opening displacement 
for selected numbers of broken fibers in panels 
having a 6061-T6 aluminum matrix. 
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Figure 37 . Load versus crack opening displacement for 89 
broken fibers in a 4-ply panel; comparison of 
the analytical and experimental curves. 
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CHAPTER V 


CONCLUSIONS 

Results of this study indicate that the mathematical 
models presented in |2] based on the classical shear-lag 
assumption predict the displacements and stresses in a uni- 
directional fibrous composite laminate with very acceptable 
accuracy for a wide range of notch lengths, laminate widths, 
and laminate thicknesses. It is also shown that the crack 
opening displacement versus load curve predicted by the 
model is very close to the actual behavior, particularly in 
the early part of the load cycle before a large number of 
fibers have broken. The extent of the stable notch exten- 
sion, that is, the number of broken but constrained fibers, 
predicted by the model is found to accurately represent that 
observed in all laminates tested. Also, the dimensions of 
the longitudinal damage region observed in the laboratory 
are found to agreo well with those calculated by the shear- 
lag model, as were the notched strengths of the test speci- 
mens . 

It is demonstrated that the thickness (number of plies) 
has very little effect on the damage growth and fracture 
behavior of a unidirectional composite laminate. Longitudi- 
nal damage regions and crack opening displacements measured 
for single-ply, two-ply, four-ply, and eight-ply specimens 
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are found to be almost Identical for a given value of remote 
strain. The stable notch extension in the form of broken 
fibers and matrix yielding ahead of the notch is found to bs 
present in the same quantities in specimens having various 
numbers of plies. 

Finally, it is shown that by heat-treating the matrix 
to increase the yield strength actually reduces the notched 
strength of the laminate. The stronger matrix transfers 
load to the unbroken fibers at a lower load than the fully 
annealed matrix, with a substantial reduction in the amount 
of matrix yielding. The results of these tests indicate 
that to increase the capability of a laminate to tolerate a 
defect during a load cycle, the matrix must be able to yield 
and thereby absorb some energy that would otherwise serve to 
propagate the notch and cause failure of the laminate. 
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